Astronomy &: Astrophysics manuscript no. 1993 


February 2, 2008 


(DOI: will be inserted by hand later) 





Discovery of magnetic fields in central stars of planetary 

nebulae"*" 

Stefan Jordan^, Klaus Werner^, and Simon J. O'Toole^ 



in 
o 
o 

(N 



> 
O 

o 

o 
in 
o 

:^ 

a.; 
6 



>< 



^ Astronomisches Rechen-Institut, Monchhofstr. 12-14, 69120 Heidelberg, Germany 
^ Institut fiir Astronomie und Astrophysik, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 
Dr.-Remeis-Sternwarte Bamberg, Sternwartstr. 7, 96049 Bamberg, Germany 

Received 10 September 2004/ Accepted 1 November 2004 

Abstract. For the first time we have directly detected magnetic fields in central stars of planetary nebulae by means 
of spectro-polarimetry with FORSl at the VLT. In all four objects of our sample we found kilogauss magnetic 
fields, in NGC 1360 and LSS 1362 with very high significance, while in EGB 5 and Abell 36 the existence of a 
magnetic field is probable but with less certainty. This discovery supports the hypothesis that the non-spherical 
symmetry of most planetary nebulae is caused by magnetic fields in AGB stars. Our high discovery rate demands 
mechanisms to prevent full conservation of magnetic flux during the transition to white dwarfs. 
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1. Introduction 

The reason why more than 80% of the known plane- 
tary nebula e (PNe) are mostly bipolar and not spherically 
symm etric (jZuckerman fc AUer Ill986l IStanghellini et akl 
Il993(l is barely understood. A popular e xplanation is the 
interacting stellar winds model ijKwok et al. 1978ib where 
the fast {v « 1000 km/sec, mass loss « lO~^M0/yr) wind 
from the central star of a PN encounters an older slow 
(v « 10 km/sec) wind from earlier phases with heavy 
mass loss (« lO~^M0/yr). The visible PN is formed in 
the shock region between both winds; if the slow wind 
was not spherical, but densest in the equatorial plane, the 
nebula is bipolar. However, neither is this model indis- 
putable, nor is the physical mechanism for the asymmetry 
of the slow wind clear. One possibility is the presence of 
a low-mass companion star which could exert a gravita- 
tional pull on the circu mstellar envelope . Rapid rotation 
and binarity (e.g. De Marco et al. Il200'3) may also cause 
asymmetries, but the most promising explanations involve 
magnetic fields. There is, however, no agreement about the 
detailed mechanism. A review on observational and theo- 
ret ical studies of the shapi ng of planetary nebulae is given 
bv lBalick FraukI ll2002l) . 

It is possible that magnetic fields from the stellar sur- 
face are wrapped up by differential rotation so that the 



later post-AG B wind will be coUimated into two lobes 
fGarcia-Scgura Another scenario says that 

magnetic pressure at the stellar surface plays a n impor- 
tant role driving the stellar wind on the AGB ijPascoli I 
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The idea that magnetic fields are important has been 
supported by the detection of polarization in radio data 
of circumstellar envelopes of AGB stars: SiO (at a dis- 
tance of 5-10 AU from the s tar), H^O (ft! 100 AU) and 
OH (100-1000 AU) masers fKemball fc Diamond I Il997t 
ISzvmczak fc Cohen lll997; .Vlemmings et al 200^ . 

For H2O masers IVlemmings et al. I l|2002() are con- 
vinced that the Zeeman interpretation is correct and that 
the magnetic field strength at the H2O maser of the Mira 
variable U Her is about 1.5G. Depending on the topology 
of the magnetic field, the corresponding surface magnetic 
field is of the order of 100-1000 G. 

The magnetic field may be either a fossil remnant from 
the progenitor on the main sequence (e.g. Ap stars), or 
can be generated by a dynamo at the interface between 
a rapidly rotating stellar core and a more slowly rotating 
envelope. iBlackman et al. I l|200l|) argue that some rem- 
nant field anchored in the core will survive even without 
a convection zone, although the convective envelop e may 
not be removed completely. iThomas et al. ] 1I1995I) have 
shown that white dwarfs which do have thin surface con- 
vection zones can support a near-surface dynamo. Since 
the field strength in their model is higher at higher lumi- 
nosities this would particularly be true for central stars of 
PNe. 
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That some central stars must contain significant mag- 
netic fields is also obvious from the fact that at least 10- 
30% of all white dwarfs have magnetic fields between 10^ 
and 10^ Gauss. Until now no magnetic fields have ever 
directly been detected in central stars of PNe. 

We have observed a sample of four central stars of plan- 
etary nebulae with high signal-to-noise (circular) spec- 
tropolarimetry between 3500 and 5900 A with the FORSl 
spectrograph of the VLT telescope . As was already demon- 
strat ed for bright white dwarfs 1 Aznar Cuadrado et al. I 
l2004l) , the unprecedented light collecting power of the VLT 
offers the possibility to investigate the presence of mag- 
netic fields on the kG level. 

2. Observations and data reduction 

The observations were obtained in service mode between 
November 2, 2003, and January 27, 2004, with the FORSl 
spectrograph of the UTl ("Antu") telescope of the VLT, 
which is able to measure circular polarization with the 
help of a WoUaston prism and rotatable retarder plate 
mosaics in the parallel beam allowing linear and circular 
polarimetry and spectropolarimetry (jAppenzeller et al~l 
Il998|l . We used grism G600B, covering the spectral range 
3400-5900 A, and a 0.8" wide slit, leading to a spec- 
tral resolution of 4.5 A. The details of our observations 
are listed in Tabled The four selected objects are bright 
{V < 12™5) and their nebulae show clear indications for 
non-spherical symmetry. So that we obtained at least one 
good result for a star, we decided to spend four times 
as much observing time on the CPN of NGC 1360 then on 
the other three objects, for which only one observing block 
was performed. In order to reduce errors from changes in 
the sky transparency, atmospheric scintillation, and vari- 
ous instrumental effects the A/4-retarder plate was rotated 
by 90° after n exposures (where n is given in the last col- 
umn of Table^. The same number of exposures were then 
taken in this configuration. 

2.1. Data reduction 

Calibration frames (bias, flat-field and He-|-HgCd arc spec- 
tra) were taken during the day, following each nights ob- 
servations. The data were reduced in the iraf environ- 
ment using the following procedure. The bias level was 
subtracted from all frames and cosmic rays were removed. 
A nightly master flat field was then constructed from each 
night's individual fiat fields. After flat-field correction, the 
stellar spectra were extracted from each frame by sum- 
ming up all CCD rows for the ordinary and extraordinary 
(e and o) beams. Background sky light was averaged over 
10 rows (giving a total of 20 rows) on either side of the 
object spectrum and subtracted. It is important to note 
that the automatic aperture and sky selection routine in 
IRAF does not always use the user-defined values, so each 
spectrum was checked manually. 

Wavelength calibration is particularly important for 
this kind of spectropolarimetric study, and special care 



was taken to ensure its accuracy. Failure to do so 
would lead to spurious polarization signals in every line. 
Calibration was done independently for the spectra of each 
beam and each position of the retarder plate (i.e. the e and 
o beams at ±45°). 

The referee suggested that spurious signals may be 
caused by using arc spectra taken at different waveplate 
angles, probably because the spectra are rebinned differ- 
ently. To test this we examined two cases: when the dis- 
persion correction was applied, all spectra were forced to 
exactly the same scale or they were simply corrected ac- 
cording to the dispersion function only. In these two cases 
the spectra were rebinned differently. When we determine 
the magnetic field strength, however, the results are the 
same within errors. This indicates that rebinning is not 
affecting our results. We have also examined the sky spec- 
tral lines at the edge of our spectra; these lines show no 
detectable polarisation, suggesting that any polarisation 
we measure is intrinsic to the star and not due to poor 
wavelength calibration. Finally we note that while instru- 
mental polarisation dominates the Stokes V/I spectrum 
when considering only one waveplate angle, we are en- 
couraged to see the polarisation profiles at the positions 
of the Balmer and Hell lines. The wavelengths are ac- 
curate to typically ~3kms~^ or ~0.05A at H/3. This is 
much lower than the spectral resolution. 

Stokes /, or unpolarized, spectra were obtained sim- 
ply by summing all spectra taken of an object in a single 
night. The Stokes V/I spectra, describing the net circu- 
lar polarization, were created by summing the exposures 
made at the same retarder plate position angle, and then 
applying the following equation 



— = , withi?^ 

/ R+1' 



fo 

fej a=+45 



X I — I , (1) 



which is equivalent to formula 4 . 1 in t he FORS 1+2 User 
Manual ("Szei fert fc Bohnhardt I I200I . Here a indicates 
the nominal value of the position angle of the retarder- 
wave plate, and fo and /e are the fiuxes on the detec- 
tor from the e and o beams of the WoUaston prism, re- 
spectively. The resulting high-quality spectra are shown 
in Fig.^ while Figs.|21-|Sl show the circular polarization 
{V/I) spectra. 

3. Determination of the magnetic field strenghs 

For weak magnetic fields (i.e. below 10 kG) theoretical 
polarization spectra {V/I) ca n be obtained by using the 
weak-field approximation (e.g.,|A ngcl fc Landstrcct " l970l: 
iLandi degl'Innocenti fc Landi degl'Innocenti ..19731) : 



V 



J -gesCzX^j-Q^{Bz) 



(2) 



where goff is the effective Lande factor (which is unity 
for Balmer lines and for the hyd rogenic He II lines; 
ICasini fc Landi degl'Innocenti lll'994l) , A is the wavelength 
expressed in A, {Bz) is the mean longitudinal component 



Stefan Jordan et al.: Discovery of magnetic fields in central stars of planetary nebulae 
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Target 


Alias 


Q 


b 


V 

(mag) 


HJD 

(+2 452 900) 


(s) 


n 


NGC 1360 


CD-26 1340 


03 33 14.7 


-25 52 18 




46.782 


104 


6 












88.725 


104 


6 












89.549 


104 


6 












90.571 


104 


6 


EGB5 


PN G211.9+22.6 


08 11 12.8 


+10 57 19 


12?^5 


88.854 


331 


3 


LSS 1362 


PN G273.6+06.1 


09 52 44.5 


-46 16 51 


12':"5 


89.796 


331 


3 


Abell 36 


PN G318. 4+41.1 


13 40 41.4 


-19 52 55 


ll'?^5 


131.773 


150 


5 



4000 7 




- 1 000 L ^ ^ ^ ^ ^ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ ^ ^ \ \ ^ ^ ^ ^ ^ ^ ^ ^ 

1000 2000 .5C0C 

input B/G 



Fig. 4. Result of the fits with simulated data having the 
same noise level as the observations of NGC 1360 for input 
magnetic fields between and 3000 G, in steps of 500 G. 
From center line to outside: mean fit result, Itr error range, 
99% confidence level, and smallest and largest fit result 
during the 1000 simulations which were performed for each 
predescribed magnetic field. 



of the magnetic field expressed in Gauss and the constant 
= e/(47rmec2) (~ 4.67 x lO'^^ Q-^ A"^). Since we do 
not have any information about the detailed field geome- 
try we can only measure the mean longitudinal field over 
the stellar surface. The maximum field strength can be 
larger than this value. 

Since both the hydrogen lines and the He II lines have 
an effective Lande factor of unity we do not expect blend- 
ing to have a large influence. However, it is clear that the 
total effect of two separate spectral lines on the polariza- 
tion is not the same as treating a blended line in the same 
way. With our method it is not possible to disentangle 
both effects. Test calculations using theoretical spectra for 
NGC 1360 have shown that the result when using blended 
Balmer and He II lines instead of a sum of non-blended 
lines (by switching hydrogen or helium, respectively, in 
the calculation of the theoretical spectrum) differ by only 
about 200 G. 



The longitudinal component of the magnetic field 
for each measurement was determined by comparing the 
observed circular polarization for an interval of ±20 A 
around the four strongest absorption lines H/?+HeII, 
Hell 4686, H7+HeII, HJ+He ll with t he prediction of 
EquationEI As in .Aznar Cuadrado et al. ( 2004) we de- 
termined (i? ^ ) by a x^-minimization procedure. Following 
IPress et al. 1 l)l986l) we determined the statistical error 
from the rms deviation of the observed circular polar- 
ization from the best-fit model. The Icr (68.3%) confi- 
dence range for a degree of freedom of 1 is the interval 
of i?2 where the deviation from the minimum is A^^ = 1; 
the 99% confidence interval corresponds to Ax^ = 6.63. 
This statistical error does not take into account any sys- 
tematic errors, particularly the blending of Balmer lines 
with Hell lines mentioned above. Only the Hell 4686 
line is not effected by blending. Although not blended, 
the weaker He II lines do not give any significant informa- 
tion; they have large statistical errors and therefore a very 
low weight. 

For each of the observation blocks, Table|21 summarizes 
our fit results for all four spectral lines and the weighted 
means Bz — Bz,iWi)/ '^Wi with i corresponding to 
the lines and Wi — 1/cr^. The total probable error is given 
by o- = (X;wj)"^^^. We list both the total Ax^ = 1 and 
Ax^ = 6.63 error range. From our statistic a significant 
magnetic field was found in three of the four NGC 1360 
observations and in the (single) observations of EGB5, 
LSS 1362, and Abell 36. However, in the latter case the 
value of the best fit is just outside the 99% confidence 
range. 

NGC 1360 clearly shows the effect of rotation between 
the observations: -1343,1708,2832, and 194 G. The dif- 
ference in time b etween the thr e e obs ervations was 42, 
0.8, and 1.0 days. IWerner et al. I l)2003j) have derived an 
upper limit for the rotational velocity of 20 km/sec from 
the with of iron lines, leading to a period larger than 0.75 
days for a radius of 0.3 R©, which is compatible with our 
result. The suc cessors of CPNs (wh ite dwarfs) are also 
rotating slowly ijKoester et al. lll998() . so that we do not 
expect any smearing out of the polarization signal during 
the observing blocks. 
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4500 



5000 
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Fig. 1. Normalized spectra of our sample of central stars of planetary nebulae (from above displaced vertically: 
NGC1360, EGB5, LSS 1362, Abcll36) 



3.1. Statistical significance of our measurements 

Since the amplitude of a polarization signal for _B « 1 kG 
is usually smaller than the la noise level of the observed 
polarization spectra, doubts about the significance of our 
result clearly originate when visually looking at the fitted 
polarization spectra. For this reason we have started a 
simulation using synthetic polarization spectra to which 
Gaussian noise of the same level as in our observation was 
added. 

In the case of NGC 1360 the noise in the single ob- 
served polarization spectra has CTnoisc = 0.0005. For given 
magnetic fields of 0, ... , 3000 G in steps of 500 G, we cal- 
culated 1000 artificial polarization "measurements" and 
treated them in the same way as our real observations. 

Figure 01 shows that the averaged weighted mean for 
the four strong spectral lines is very close to the given 
value of the magnetic field. It also shows that for an as- 
sumed magnetic field of B = G, only one result reaches 
900 G. If we conservatively assume that the systematic er- 
ror is 500 G, two of the four observing blocks of NGC 1360 
have a much larger measurement (1708 and 2822 G), and 
one (-1342 G) is only marginally below this extremely pes- 
simistic criterion. Of all the simulations, 99% with an as- 
sumed B = OG have fitted field strengths below 660 G. On 



the other hand, if we assume a magnetic field of 1000 G, 
the fits to the artificial spectra result in values between 
280 G and 2110 G, with 99% of them lying between 342 
and 1670 G. 

The lower panels in Fig|3 show an example of one 
of the 1000 artificial spectra for an assumed magnetic 
field of 1500 G, which is close to the 1708 G value for 
NGC 1360 measured from the second observing block. 
It makes it clear tha t , as already demonstrated by 
^.znar Cu adrado et al. I l(2004l) . visual inspection is mis- 
leading, since the eye does not take into account an av- 
erage small excess of right- and left-handed polarization 
on different sides of the line core, respectively, which con- 
tributes to our analysis. The standard deviation of all 
1000 fits is 254 G, very close to our formal la error for 
1708 G, which is 258 G. We therefore conclude that in the 
case of NGC 1360 the statistical errors from our anal- 
ysis are indeed realistic in order to judge how accurately 
the magnetic field can be determined. 

A somewhat different situation occurs in the case of 
LSS 1362 {Bfit — 1891 G), where the noise level of Cnoise = 
0.00085 is larger. If we assume that no magnetic field ex- 
ists, four of the 1000 simulations result in a fitted magnetic 
field strength exceeding 1891 G. If, in order to account for 
a possible systematic error, we set the limit at 1500 G, 16 
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Fig. 2. Circular polarization [V/I) in the four observation blocks of of the central star of NGC 1360 in the vicinity 
of the strong spectral lines H(5+HeII, H7+HeII, Hell 4686, H/3+HeII compared to the prediction by the low-field 
approximation (Eq.|21) using a longitudinal magnetic field of — 1343G, 1708 G, 2832 G, and 194 G, respectively 

Table 2. Magnetic fields derived from the four strongest lines in our sample of central stars of planetary nebulae. The 
error margins correspond to a la (68.3% confidence) and 6.6cr (99% confidence) level. 



Target 


Date 




B/G 




B/G 






H(5+HeII 


H7±HeII 


Hell 4686 


H/3+HeII 


total 


NGC 1360 


03/11/03 


-493 ± 835 


2483 ± 688 


-1114 ±427 


-1355 ±413 


-1343 ±259 






(±2153) 


(±1772) 


(±1101) 


(±1066) 


(±668) 


NGC 1360 


14/12/03 


342 ± 624 


3553 ± 528 


1726 ± 426 


768 ± 543 


1708 ± 257 






(±1607) 


(±1361) 


(±1099) 


(±1399) 


(±664) 


NGC 1360 


15/12/03 


3146 ± 735 


2082 ± 695 


3714 ± 389 


1324 ± 557 


2832 ± 269 






(±1895) 


(±1790) 


(±1002) 


(±1437) 


(±695) 


NGC 1360 


16/12/03 


2548 ± 1024 


1303 ± 505 


-1176 ±591 


316 ±438 


194 ± 277 






(±2638) 


(±1302) 


(±1522) 


(±1128) 


(±2548) 


EGB5 


14/12/03 


577 ± 1171 


2875 ± 1022 


484 ± 4707 


2171 ±837 


1992 ± 562 






(±3016) 


(±2633) 


(±12121) 


(±2156) 


(±1449) 


LSS 1362 


15/12/03 


-299 ± 1295 


2089 ± 729 


3608 ± 894 


1550 ± 531 


1891 ± 371 






(±3335) 


(±1878) 


(±1848) 


(±1368) 


(±912) 


Abell 36 


26/01/04 


1863 ± 1144 


-448 ± 949 


1842 ± 1134 


1553 ±719 


1169 ±466 






(±2946) 


(±2444) 


(±2920) 


(±1852) 


(±1202) 
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Fig. 3. Circular polarization {V/I) observed in the central stars of LSS 1362, Abell36, and EGB 5 in the vicinity 
of the strong spectral lines H5+HeII, H7+HeII, Hell 4686, H/3+HeII compared to the prediction by the low-field 
approximation (Eq.EJ using a longitudinal magnetic field of 1891 G, 1169 G, and 1992 G. In the lower level we show an 
example for a fit to one of the artificial spectra with a noise level of 0.0005 and an assumed magnetic field of 1500 G. 
The fit results in B = 1763G 



(1.6%) of the simulations provide a larger field strength. 
The standard deviation for an assumed field strength of 
2000 G is 646 G, about 75% larger than the formal Icr er- 
ror from the analysis. Therefore the probability that 
LSS 1362 has a magnetic field of more than 1000 G is very 
high. 

In the case of Abell 36, where we measured a magnetic 
field of 1169 ± 466 G, the situation is more uncertain: for 
Cnoise = 0.00067 wc find that for an assumed magnetic 
field of G 144 (14.4%) of all artificial polarization spectra 
mimic a magnetic field > 1169 G, 555 (55.5%) a magnetic 
field larger than 669 G (if we again estimate the maximum 
systematic error to be 500 G). Therefore, we would not 
regard the derived magnetic field as very significant. 

Although we formally measured a magnetic field of 
1992 ± 562 in EGB 5, the case for a kilogauss magnetic 
field is probable but not with the high certainty indi- 
cated by the error range from the analysis. For G and 
fnoise = 0.0012 we find that 64 (6.4%) models exceeded 



1992 G, and 142 (14.2%) the limit of 1492 G, taking into 
account systematic uncertainties. Due to the higher mea- 
sured value, this is a clearer case than that of Abell 36. 

Our simulations with artificial polarization spectra 
clearly show that much more realistic error estimations 
can be obtained compared to the formal errors from the 

analysis. They show, however, that our determinations 
of magnetic fields are significant in the case of NGC 1360 
and LSS 1362 even though the maximum polarization sig- 
nal does not exceed the noise level. 

4. Parameters of the target stars and nebulae 

Atmospheric parameters for the central stars of 
NGC 1360, Aben36, and LSS 1362 were derived by 
iTraulsen et al. I ^l2005^ from HST STIS and optical spec- 
tra. For EGB 5 T pfr and logg were determined by 
iLisker et al. I l|2004j) . From these values, masses and radii 
for the central stars were estimated, as well as the masses 
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and radii on the main sequence and of the white dw arf suc- 
cessor s. For this purpose, mass-radius relations bv lWoodI 
were used. In the case of EGB 5 no such values 
could be derived, since its ce ntral star (a hot s ubdwarf ) is 
a result of binary evolution l|Karl et al. 112003^ . The value 
of these parameters together with a designation of the 
planetary nebular morphology is listed in Table|31 

If we assume complete conservation of magnetic flux 
through the stellar surface from the main sequence to 
the white dwarf stage, we can estimate the magnetic field 
strength of the precursors and successors. The magnetic 
field strength measured from the third observing block in 
the central star of NGC 1360 of 2800 G would translate 
into a field strength of 50 G on the main sequence while 
the field strength will be enhanced to 2 MG is the star will 
reach the white dwarf stage. For Abell36 (1170 G) and 
LSS1362 (1900 G) the values are 9.3 G, 0.35 MG, 24 G, 
and 0.43 MG, respectively. 

This is surprising, because magnetic fields of 0.35- 
2.0 MG would be detectable from Zeeman splitting in 
high-resolution an d high-signal-to-noi s e spec tra, e.g. from 
the SPY survey ijNapiwotzki et al. 1 12003|) and in the 
majority of the sample stars such high magnetic field 
strengths can be excluded. Therefore, we have to assume 
that our assumption of full conservation of magnetic flux 
is invalid. This might be a hint that the magnetic field is 
not strongly concentrated to the degenerate stellar core, 
where the t ime scale for the decay should be of the order of 
10^° years (jChanmugam fc Gabriel 1ll972HFontah^ et al. I 
Il973l) . It could instead be present in the envelope, where 
it might be destroyed by convection or mass-loss. 

5. Discussion and conclusions 

We have detected magnetic fields in 50%-100% of our 
small survey for magnetic fields in central stars of plan- 
etary nebulae, depending on how conservatively the cri- 
teria for statistical significance are set. This provides 
very strong support for theories which explain the non- 
spherical symmetry (bipolarity) of the majority of plan- 
etary nebulae by magnetic fields. In this first survey we 
have not performed a cross check with any spherically- 
symmetric nebulae, although this is planned as a follow- 
up. 

Although based on only four objects, our extremely 
high discovery rate demands that magnetic flux must be 
lost during the transition phase between central stars and 
white dwarfs: if the magnetic flux was fully conserved, 
our four central stars will have fields between 0.35 and 
2 MG when they become white dwarfs. Although the num- 
ber of white dwarfs with magnetic fields is still a mat- 
ter of debate, with a range between about 3 and 30%, 
even the latter value, which includes object s with kG 
field strengths fAznar Cuadrado et al. l l2004j) . is far off 
our high number. Liebcrt et al. . (2003^1 quantified the in- 
cidence of magnetism at the level of ~ 2 MG or greater to 
be of the order of '^10%. This argument would not change 
by much if we consider that we have so far only looked at 



central stars with non-spherical symmetric nebulae. An 
almost 100% probability of magnetic fields larger that 
100 kG can be excluded by t he data from the SPY sur- 
vev llNapiwotzki et al. l2003^ as well as the sample from 
lAznar Giiadrado et al. l)2004l) . It is also worth mention- 
ing that our central stars have typical white dwarf masses 
(0.48-0.65 Mq) and are not particularly massive. White 
dwarfs with MG f ields tend to be more massive than non- 
magnetic objects l|Liebert Ill988l) . 

If the magnetic field is located deep in the degener- 
ate core of the central star, it is very difficult to imag- 
ine a mechanism to destroy the ordered magnetic fields. 
Therefore, it would be more plausible to argue that the 
magnetic field in the central stars is present mostly in 
the envelope where it can be affected by convection and 
mass-loss. For central stars hotter than 100 000 K we do, 
however, not expect convection; only in the central star of 
EGB 5 we cannot exclude such a mechanism. 

If we assume that the magnetic fields are fossil and 
magnetic flux was conserved until the central-star phase, 
we estimate that the field strengths on the main sequence 
were 9-50 G, which are not directly detectable. Therefore, 
our measurement may indirectly provide evidence for such 
low magnetic fields on the main sequence. 

Polarimetry with the VLT has led to discovery of mag- 
netic fields in a large number of objects in the final stage 
of ste llar evolution: white dwarfs l| Aznar Cuadra do et al. I 
hot subdwarf stars l|0'Toole et al. .^2005) . and now 
in central stars of planetary nebulae. Although we have 
now provided a good basis for the theoretical explanation 
of the planetary nebula morphology - which can more 
quantitatively be correlated with additional observations 
in the future - new questions about the number statistics 
of magnetic fields in the late stages of stellar evolution 
have been raised. 
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